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The conductivity of doped graphene is considered taking into account scattering by short-range 
nonsymmetric defects, when the longitudinal and transverse components of conductivity tensor 
appear to be different. The calculations of the anisotropic conductivity tensor are based on the 
quasiclassical kinetic equation for the case of monopolar transport at low temperatures. The effective 
longitudinal conductivity and the transverse voltage, which are controlled by orientation of sample 
and by gate voltage (i.e. doping level), are presented. 

PACS numbers: 72.10.Fk, 72. 20. Dp, 72.80.Vp 



m 



CD 



cd 



O 

o 



> 

00 

en 

in 

o 
o 



X 



The conductivity tensor aap is determined both by the 
symmetry properties of carriers and by the character of 
scattering processes. For a medium with cubic (quadratic 
for the 2D case) symmetry, the conductivity appears to 
be a scalar, i. e. aaj3 oc SajS- [ij Anisotropic conductivity 
of hot electrons in low-symmetric bulk materials (Ge and 
Si) were studied over 50 years ago. '2| For the hexago- 
nal symmetry case, which is correspondent to an ideal 
graphene sheet, the longitudinal (along X-axis, see Fig. 
1) and transverse conductivities of hot carriers appear to 
be different. This anisotropic conductivity appears due 
to the trigonal warping of energy spectrum [3| even for 
the isotropic scattering case. Such an anisotropy is weak 
[oc {ap/%Y, where a is the lattice constant and p is the 
characteristic momentum] and it can be essential for the 
high-energy carriers, at energies > 1.5 eV. In the linear 
regime, the anisotropy of conductivity is possible due to 
scattering by short-range nonsymmetric defects, see Fig. 
lb where a substitutional impurity is shown. Although 
the nonsymmetric defects were discussed in Refs. 4-7, 
the anisotropy of conductivity was not theoretically an- 
alyzed or mentioned. To date, no experimental data are 
reported concerning any mechanism of anisotropy listed; 
so that, a study of this phenomena is timely now. 
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FIG. 1: (Color online) (a) Equal enegry contours around the 
K and K' points with step 0.5 eV; a visible anisotropy takes 
place at energies > 1.5 eV. (b) Graphene sheet with a short- 
range nonsymmetric defect marked by ellipse. 

In this Letter, we consider the low-temperature 
conductivity taking into account both the long-range 
isotropic disorder (see |8|4lO|) and the short-range non- 
symmetric defects (see [3t31)- The longitudinal and 



transverse components of conductivity tensor are calcu- 
lated based on the linearized Boltzmann equation with 
the transition probabilities written in the Born approx- 
imation. The effective conductivity and the transverse 
voltage of a graphene strip are presented for the case of 
weak anisotropy and their dependencies on carrier's con- 
centration (gate voltage) and on orientation of strip are 
discussed. 

We consider the momentum relaxation in graphene 
which is caused by the interaction of carriers with a weak 
potential ^ ■ ?7(x — x^), where Xj is the coordinate of 
jth impurity (j = 1, . . .Nim and Nim is the number of 
impurities over the normalization area L^). The poten- 
tial of a single impurity placed at x = is given by 
?7(x) = Ux + u(5(x); here we separated the long-range 
scalar contribution, Wx, and the short-range addend de- 
termined by the 4x4 matrix u. For a nonmagnetic de- 
fect, this matrix should be Hermitian and time- reversal 
symmetric and one can write u using 9 independent real 
parameters Ugi: 
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where 4x4 matrices Sg and A/ are introduced in Refs. 
4-6. In the framework of the Born approximation, the 
transition probability between the c-band states jjy'p') 
and jryp) is given by 
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where p is 2D momentum, nim = Ni,-n/ L"^ is the impurity 
concentration, rj — K,K' is correspondent to K- or K'- 
valleys and the linear dispersion law Sp = v\yP is not 
dependent on 77 {vw = 10^ cm/s is the characteristic 
velocity). ^^ Below we restrict ourselves by the case of 
weak short-range corrections to the long-range scattering 
contribution. As a result, one obtains the intravalley 
matrix element in (2) as 
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where the main contribution is written through the 
Fourier component of Wxi which is dependent on the mo- 
mentum transfer Tiq = p — p', and through the factor 
(1 — cosA(^) with A(p — p',p described the suppression 
of backscattering processes. The A(^-dependent contri- 
butions, which are proportional to u^;, were omitted in 
Eq. (3). In addition, a weak A(^-dependent intervalley 
matrix element, which is proportional to u^; or to w^j 
with (s, I) —i' {x, y) should also be neglected. The re- 
maining corrections in Eq. (3), which are dependent on 
(/?-!-(/?' {ip and ip' are the polar angles of p and p'), are 
responsible for the anisotropic conductivity under con- 
sideration if M^2 7^ M^2i i-6- scattering processes along 
X- and y-directions are different. 

Next, we solve the Boltzmann kinetic equation lin- 
earized with respect to the steady-state electric field E. 
Taking into account that the transition probability given 
by Eqs. (2), (3) does not dependent on rj and that 
Wpp' — M^p'p, we write the kinetic equation for the 
anisotropic part of distribution A/p, which is ex E, in 
the following form: 



e(E-Vp)^ 
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Here Vp = VpSp is the velocity, /^ is the equilibrium 
distribution and dfe/de ~ —5{£f —£p) for the degenerate 
carriers with the Fermi energy ep. We solve the integral 
equation (4) using the variational approach |ll| with the 
trial distribution function A/p — e (C • Vp) S{eF — £p). 
An unknown vector C is determined from the extremum 
conditions for the quadratic form: 



/C(C) = (C-i-E +(C-L-C), 



(5) 



where the field-dependent contribution is determined by 
the tensor: 



Aap = -yf Yl '^('^-P " ^P) 



—=0al3 -. (Dj 
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written through the density of states at the Fermi en- 
ergy, pp. The collision-induced contribution in Eq. (5) 
is determined by the matrix: 



part of potential (vq = Wg=o) with the correlation lenth Ic- 
The function "^(pplc/h) was introduced in Q for the case 
of the finite-range potential with Gaussian correlations. 
The current density is given by I = 
(4e/L^)^ VpA/p, where the factor 4 is due to 
the spin and valley degeneracy. Since A/p oc C oc E, 
the conductivity tensor a is introduced as I = itE. 
Substituting A/p into I and performing integration 
over p, one obtains I — [evwYCpp /2, where C is 
determined from the extremum conditions which give 
the linear equations: dlC{C)/dCa = for a = x,y. As a 
result, one obtains the diagonal tensor aaa — o" ± S(t/2, 
where -|- (or — ) is correspondent to axx (or ayy). Here 
the isotropic part of conductivity a and the anisotropic 
correction 6a are given by 
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so that the anisotropy of conductivity is determined by 
the dimensionless factor {u^^ — Uy^)/(8vQ). 
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FIG. 2: (Color online) (a) Graphene strip of width d with an- 
gle 6 to X-direction (dotted arrows are correspondent to the 
XQY coordinate system). Longitudinal and transverse fields, 
E± and _E||, as well as longitudinal current /y are shown, 
(b) Functions ^~^ and ^"^ which describe gate voltage de- 
pendencies of transverse voltage, E±/E\\, and of anisotropic 
contribution to conductivity, {o^ff — a), given by Eq. (10). 

Using the conductivity tensor determined by Eqs. (9), 
we consider the conductivity of strip, see Fig. 2a where 
d is the width of strip and 9 is the angle between its ori- 
entation and X-axis. Since the absence of the transverse 
current, I± = 0, the linear relation I = ctE permits one to 
calculate the effective conductivity, which is introduced 
where -I- (or — ) is correspondent to L^x (or Lyy). The according to In = aeffEu, and the induced transverse 
factor Vq'${ppIc/K) appears here due to the long-range voltage, E±d. For the weak anisotropy case, we obtain 



and after substitution of the matrix elements (3) one ob- 
tains Lap = if a 7^ /?. The diagonal components of Eq. 
(7) take form: 
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the simple angle dependencies 
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■ sin 26*, (Te// — cr 
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and maximal values of the transverse voltage |i<^x/£^||| 
and of the anisotropic contributions to (Je// appear at 
— 7r/4, 37r/4, . . . and at 6* = 0, 7r/2, . . ., respectively. 
The concentration (or gate voltage, Vg) dependencies are 
determined through the function "^{pplc/K) introduced 
in Eqs. (8), (9). According to Eq. (10), one obtains 
|i?j_/i?ll| oc ^^^ and 5a ex \&~^; we plot these functions 
in Fig. 2b for the graphene strip placed over the Si02 
substrate of 300 nm width. We also use Ic —10 nm in 
agreement with experimental data, see Ref. 9. Since the 
function ^ decreases with V^, both the transverse voltage 
|i?j^/-B||| and the anisotropy of conductivity {(Jeff — ct) 
increase with Vg] notice, that the anisotropy of conduc- 
tivity increases more faster (about 10 times at Vg ~10 
V) because Sa ex \I>~^. The transverse field and the 
anisotropy of effective conductivity given by Eq. (10) are 
proportional to the unknown parameter (u^^ — u^^)/(8wq) 
and a value of the anisotropy effect is not estimated here. 
However, these results permit one to verify a mechanism 
of anisotropy using the gate voltage and angle dependen- 
cies obtained. In addition, a ratio [u^^ ~ '^^z) I {^"^V) '^^'^ 
be measured. Supposing [u^^ — u'^^)/{Sv'^) ^ 10~^, one 
obtains that \E±/E\\\ and \Sa/a\ are about 10 %. 

Further, we discuss the assumptions used. The con- 
sideration performed is based on the simplified descrip- 
tion of scattering with separated long- and short-range 
parts of potential. The parameters of long-range scatter- 



ing are determined from the phenomenological consider- 
ation [9J while the short-range contribution is expressed 
through unknown matrix it, without verification of a mi- 
croscopical nature of defects. The variational solution of 
Eq. (4) gives a good estimation of the anisotropy un- 
der consideration. The results (10) do not dependent 
on a width of strip if d exceeds the mean free path due 
to the long-range scattering, so that an edge scattering 
is not essential under the condition d 3>0.1 /xm (here 
we used the results for non-ideal nanoribbons [125). We 
also restrict ourselves by the degenerate carriers case be- 
cause the anisotropy effects increase with concentration 
(or gate voltage). The rest of assumptions (quasiclas- 
sical kinetic equation, weak acoustic phonon scattering, 
and model description of long-range scattering [9|) are 
rather standard. 

In conclusion, we have demonstrated that the nonsym- 
meric short-range scattering gives rise to the transverse 
voltage and to the orientation-dependent effective con- 
ductivity of graphene strip. These dependencies permit 
one to measure asymmetry of scattering which is ex- 
pressed through the characteristics of short-range defect. 
We believe that these results will stimulate further mea- 
surements of anisotropy and microscopical calculations 
of nonsymmetric defects. Beside this, a nonequivalent 
heating of carriers in different valleys (and an intervalley 
redistribution of carriers) by a strong electric field takes 
place under an essential anisotropy of scattering (in anal- 
ogy with the bulk case [13]). 

I am grateful to E. I. Karp for a useful comment. 
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